home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libconv / TRY / simple_diir2d.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  702 b   |  34 lines
  1.     subroutine simple_diir2d( nfx,nfy, ngx, ngy, f, g, h)
  2.     double precision f(0:(nfx-1),0:(nfy-1))
  3.     double precision g(0:(ngx-1),0:(ngy-1))
  4.     double precision h(0:(nfx-1),0:(nfy-1))
  5.     integer nfx, nfy, ngx, ngy
  6.  
  7.     integer i, j, ih, jh, i1, i2, j1, j2
  8.  
  9.     do j = 0, nfy-1
  10.         do i = 0, nfx-1
  11.         h(i,j) = f(i,j)
  12.         end do
  13.     end do
  14.  
  15.     do jh = 0, nfy-1
  16.         do ih = 0, nfx-1
  17.         j1 = max( 0, jh - (ngy-1))
  18.         j2 = jh - 1
  19.         i1 = max( 0, ih - (ngx-1))
  20.         i2 = ih - 1
  21.         do j = j1, j2
  22.             do i = i1, ih
  23.             h(ih,jh) = h(ih,jh) - h(i,j) * g(ih-i, jh-j)
  24.             end do
  25.         end do
  26.             do i = i1, i2
  27.             h(ih,jh) = h(ih,jh) - h(i,jh) * g(ih-i, 0)
  28.             end do
  29.             h(ih,jh) = h(ih,jh) / g(0,0)
  30.         end do
  31.     end do
  32.     return
  33.     end
  34.